load my_signal_with_noise5.mat;
% 傅里叶变换  
Y = fft(signal_noisy);
f = (0:length(Y)-1)*(fs/length(Y));

% 频率掩膜，去除高频噪声
mask = (f < 100); % 假设信号主频率在100Hz以下
Y_filtered = Y .* mask;

% 逆傅里叶变换
y_denoised = ifft(Y_filtered, 'symmetric');

% 绘图
figure;
subplot(3,1,1);
plot(t, signal,'b');
hold on;
plot(t, y_denoised,'r');
hold off;
title('蓝色为原始信号，红色为处理后的信号');
xlabel('时间 (s)'); ylabel('幅度');

subplot(3,1,2);
plot(t, signal_noisy);
title('带噪声信号');
xlabel('时间 (s)'); ylabel('幅度');

subplot(3,1,3);
plot(t, y_denoised, 'r');
title('频率掩膜后的信号');
xlabel('时间 (s)'); ylabel('幅度');
